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2.3 A gravity current 



For the highly nonhnear equation, a relatively simple solution is that of a stationary (or 
permanent) wave which is profile advancing at a constant speed without changing its shape. 
Mathematically the profile is describable as 



/i(a;, t) = h{x — Ct) = h{a), a = x — Ct 
By the chain rule of differentiation, 

dh{x — Ct) dh da dh 



(2.3.1) 



dh{x — Ct) dh da ^dh 
dt da dt da'' 



dx 



da dx da 



Hence (2.3.30) reduces to an ordinary differential equation. 
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(2.3.2) 



Integrating once we get 
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Let the gravity current advance along a dry bed, then /i = is a part of the solution. The 
constant of integration must be set to zero. Introducing the dimensionless variables 



h = Hch', a — Lea' , with Lc = He/ tan9, 
where He is the maximum depth far upstream, we get 
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(2.3.3) 



(2.3.4) 



Let the gravity current be uniform far upstream, then 

dh' 



h' 1, 



da' 



0, as a' — >• — oo. 



(2.3.5) 
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It follows that 



= 1 



pgH^ sin 9 



or, 



pgii?^ sin ^ 



(2.3.6) 
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and 



(2.3.7) 



One of the solution is h' — 0, representing the dry bed. For the nontrivial solution, we 
rewrite 



This is an implicit relation between h' and a', and represents a smooth surface decreasing 
monotonically from h = 1 at a' ^ — oo to h' = at the front a' = a'^, as plotted in Figure 
2.3.1. Note from (2.3.8) that da'/dh' = when h' = 0, implying infinite slope at the tip of 
the gravity current. This infinity violates the original approximation that dh'/da' = 0(1). 
Fortunately it is highly localized and does not affect the vahdity of the theory elsewhere (see 
Liu & Mei, 1989, JFM). 




(2.3.8) 



which can be integrated to give 




(2.3.9) 



Figure 2.3.1: Gravity current down an inclined plane 



Eq. (2.3.6) tells us that the speed of the front is higher for a thicker layer, steeper slope 
or smaller viscosity. This relation can be confirmed by a quicker argument. In the fixed 
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frame of reference, the total flux must be equal to CH. therefore C must be equal to the 
depth-averaged velocity u which is given by (2.3.19) with dh/dx — 0. 

A similar analysis has been applied to a fluid-mud which is non-Newtionian characterized 
by the yield stress. Laboratory simulations have been reported by Liu & Mei (J. Fluid Mech. 
207, 505-529.) who used a kaolinite/water mixture. Figure 2.3.2 shows the setup of the 
inclined flume and Figure 2.3.3 shows the recorded profiles of the gravity current along with 
the theory . The agreement is very good, despite the steep front where the approximation 
is locally invalid. 



Figure 2.3.2: Experiment setup for a mud current down an inclined plane. From Liu & Mei 
1989. 
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Figure 2.3.3: Profiles of a mud current down an inclined plane. Prom Liu & Mei 1989. 



